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Abstract 

We discuss in a systematic way the extraction of neutrino masses, mixing angles and leptonic 
CP violation at neutrino factories. Compared to previous studies we put a special emphasis 
on improved statistical methods and on the multidimensional nature of the combined fits 
of the v e — > u^, i> e — > appearance and — > v^, — + disappearance channels. Uncer- 
tainties of all involved parameters and statistical errors are included. We find previously 
ignored correlations in the multidimensional parameter space, leading to modifications in 
the physics reach, which amount in some cases to one order of magnitude. Including proper 
statistical errors we determine for all parameters the improved sensitivity limits for various 
baselines, beam energies, neutrino fluxes and detector masses. Our results allow a compar- 
ison of the physics potential for different choices of baseline and beam energy with regard 
to all involved parameters. In addition we discuss in more detail the problem of parameter 
degeneracies in measurements of 5cp- 
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1 Introduction 



The potential to measure neutrino masses, mixings, matter effects and leptonic CP violation 
at neutrino factories has already been studied for different setups (see e.g. [1, 2, 3]). Among 
the main points discussed are the achievable precision for the oscillation parameters, the 
optimal baseline and the best muon energy, the search for matter effects, the search for CP 
violation and specifically, whether measurements of the CP phase are possible. The results 
of such studies depend, however, strongly on the chosen experimental setup and the assumed 
physics parameters, and the combined total parameter space is not easily overlooked. 

A common strategy to deal with this situation is to discuss only one or two parameters 
simultaneously, while the remaining parameters are set to "standard values" (like the muon 
flux or sin 2 2^23 = 1)- This is a good method for first estimates of the physics capabilities 
of such experiments. Statistical errors computed in this way are, however, underestimated 
since they miss possible correlations with the parameters which are hold fixed. This can 
lead to quite severe errors, especially close to the sensitivity limits, where a strong interplay 
between several parameters exists. Similar problems arise in a two step analysis, where 
the disappearance channel is analyzed with sub-leading corrections ignored, and where the 
extracted leading parameters are used for the extraction of sub-leading parameters. 

We present in this work an event rate analysis based on a statistical method, which treats all 
parameters on equal footing instead of arbitrarily selecting specific parameter subsets. The 
best method would be to perform full six parameter fits to the simulated event rates, but this 
is not feasible, since the required computing time for the exploration of the parameter space 
and the parameter dependencies is enormous. We adopt therefore a somewhat simplified, 
but still general method which is based on the calculation of all two dimensional slices 
through the fit manifold. This method is related to the covariance matrix method and it 
automatically includes all two parameter correlations. We apply this method to the analysis 
of all involved physics parameters as a function of the machine parameters. In particular, 
we give improved sensitivity limits for measurements of the mixing angle #i 3 , matter effects 
and CP violation. There, the impact of the improved statistical treatment is most evident: 
Sensitivity limits can be deteriorated by up to one order of magnitude by correlations with 
sub-leading parameters. We demonstrate furthermore how input from other experiments 
(like KamLAND [4]) helps to improve the results obtained from neutrino factories. Finally, 
we study the possibility of combining two distinct baselines and we present results, which 
help to determine the optimal baseline and the best muon energy. 

The paper is organized as follows: First we present in section 2.1 analytic formulae for 
the relevant oscillation probabilities in vacuum and in matter, which allows to understand 
many numerical results qualitatively. In section 2.2 we provide a classification scheme of 
the parameter space, which is important, since it provides an overview of how we analyze 
and discuss the complex parameter space in detail in section 4. The framework which we 
use for our numerical simulation of event rates has been discussed in detail in earlier works. 
We explain this framework therefore only briefly in section 2.3 and we give references to 
earlier studies. The statistical methods which we used in our study are described in detail 
in section 3. The results are given in section 4, where we follow the classification scheme 
introduced in sec. 2.2. Finally, in section 5 we conclude. 
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2 Framework 



2.1 Oscillation probabilities in vacuum and in matter 

We assume standard three neutrino mixing with the leptonic mixing matrix U, which coin- 
cides with the standard parameterization of the quark mixing matrix [5] . The probabilities 
which describe in vacuum flavour oscillations between an arbitrary number of neutrinos are 

Pi m = S lm -4j2 M4h ei ) sm 2 A kh -2J2 M^D sin(2A fcft ) , (1) 

k>h k>h 

where Jg ej := U^U^.U^^, A kh := Am^/(4^) and A := A 31 . Matter effects 
lead to sizable corrections of these probabilities. The numerical results shown in section 4 are 
based on event rates, calculated from the full oscillation formulae, with the matter potential 
included. As usual, it is assumed that only muons can be detected and a number of the 
results can be understood analytically from the corresponding probabilities P{f^ — > Vp) and 
P{y e — > Vfj). The full expressions are, however, quite lengthy and do not allow much insight. 
We provide therefore simplified expressions by expanding these probabilities in the small 
mass hierarchy parameter a := Am^/Am^. Interesting CP- and Am 2 2 -effects will only 
occur for bi-maximal mixing and if a is not too small, i.e. for the LMA-MSW solution of the 
solar neutrino deficit. We focus therefore on the currently favoured LMA-MSW solution 
of the solar neutrino deficit [6]. The expansion in a is a good approximation as long as 
the oscillation governed by the solar Am^ is small compared to the leading atmospheric 
oscillation, i.e. «A < 1. This translates into a upper bound for the baseline: 

£ < 800 0km f^f«. (2) 



v GeV J V Am 21 ... 

Up to order a 2 one obtains for the vacuum appearance and disappearance probabilities: 

Piv e -> v^) w sin 2 6*23 sin 2 29 13 sin 2 A 

± a sin 5qp cos #13 sin 2# 12 sin 2# 13 sin 2# 23 sin 3 A 

+ a cos 5cp cos #13 sin 29 i2 sin 29i 3 sin 2# 2 3 cos A sin 2 A 

+ a 2 cos 2 # 23 sin 2 2# i2 sin 2 A (3) 

P(h>^ — > v^) ps 1 - cos 2 6*i3 sin 2 2# 23 sin 2 A 

+ 2a cos 2 0i3 cos 2 0i2 sin 2 26> 23 A cos A (4) 

The numerical magnitude of the different terms in the expansion is affected by different 
powers of the small mixing angle $13. All terms of the expansion can be written in the form 
a na 9™3 and for a ~ 6> 13 the numerical size of each term is roughly controlled by n = n a + n g . 
This amounts effectively to a reordering of the expansion in eq. (3), where all terms have 
n — 2. Note that all terms proportional to a have also a sin 2#i3 factor and are thus typically 
of the same magnitude as the a 2 -term. Higher orders of the expansion in a are, however, 
always suppressed relative to these leading terms, since rig > 0. The precise magnitude of 
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the terms which are in the sense of this reordering of the same magnitude depends of course 
on the parameter values, especially on the size of a compared to 9 13 . For 9 13 close to the 
present upper bound (sin 2 2# 13 0.1), the first term of eq. (3) is, for example, dominating 
and the last term which is proportional to a 2 is tiny and can be ignored. For smaller 
values of #13, all four terms have approximately the same importance and all of them have 
to be considered for analytical explanations. The last a 2 -term in eq. (3) can only become 
dominating for extremely tiny values of 13 . 

The CP phase <5cp produces only in the appearance channel significant effects. This can be 
seen from eq. (4), where 5qp does not show up, while 5cp is contained in eq. (3). Both terms 
in eq. (3), which contain information on the CP phase are, however, suppressed by the mass 
hierarchy. The precise value of Am^ has thus considerable impact on the magnitude of CP- 
violating effects. There is a profound difference between seeing CP-violating effects directly 
in an experiment and optimally measuring the CP phase 5qp on the other hand. The point 
is that the term in (eq. 3) which is proportional to sin Sqp changes sign when anti-neutrinos 
are considered and violates CP explicitely. The CP phase enters, however, also in the term 
which contains cos<5cp, which does not violate CP explicitely, but still can and should be 
used as an important lever arm to extract 5qp from measurements. We treat the CP phase 
5qp in our analysis therefore exactly in the same way as all other parameters, i.e. we fit 
rates to the full equations including the cos^cp-term and constrain the parameter space 
of Sqp- Our numerical analysis also includes matter effects naturally, and the procedure 
has therefore the important advantage that it maximally exploits the information on 5qp 
contained in the appearance channels 1 . 

For baselines above some hundred kilometers, one must include also the matter potential felt 
by neutrinos passing through the Earth. This opens the interesting possibility to observe 
MSW-effects [7, 8, 9, 10] in the appearance channel, while the disappearance probability 
P{ v n ~^ v \i) is on ly marginally affected by matter effects. We include these effects in our 
numerical analysis and for an analytic discussion we are interested in a reliable and traceable 
expression for P[y e — > z/ M ) in matter of constant average density. A simple approximative 
result which is valid for small values of #13 was derived in [11], where also expressions for 
larger values of #13 are given. Similar formulas were also derived in [12, 2]. 

P(v e ^v\ « sin 2 g 23 sin 2 2g 13 Sin2[(i ~ 1)A] 

V e til Ai 13 _ ^ 

sinfAA) sinlYl — A)A1 

± a sin<5cpcos#i3sin26> 12 sin26> 13 sin26>23sin(A) 



+ a cos 5qp cos #13 sin 26 12 sin 2#i 3 sin 2# 23 cos (A) 



A (1-A) 

sin(iA) sin[(l - A) A] 
~~A (1-A) 



+ a 2 cos 2 g 23 sin 2 2g 12 Sin2 j^ A) . (5) 



1 Note that seeing explicit CP violation is not easy in the presence of matter. In vacuum an asymmetry 
between neutrino and anti-neutrino appearance rates would be a clear signal for CP violation, but such 
asymmetries arise in matter also from MSW effect, such that both effects must be carefully separated. 
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Here, A = Aj Am\ x = 2VE U / Am^ and V = \p2Gpn e where Gp is the Fermi coupling 
constant and n e the electron density of the involved matter profile. The expressions show 
at first sight that in the limit of small baselines, the vacuum result (eq. 3) is recovered. 

2.2 Classification of oscillation parameters 

The full three neutrino oscillation formulae eq. (1) have in general a rather complex pa- 
rameter structure. All mass splittings and mixing parameters, namely Am^, Aro^, #12, 
#13, #23 and Sep appear in the relevant transition probabilities even after the expansion in 
a of section 2.1. The expansion allows, however, to make use of two experimentally justi- 
fied facts: First, the hierarchy of the neutrino mass splittings, i.e. a = Am^/Am^ <C 1, 
which allows to identify hierarchy suppressed "small Am 2 effects" . Second, the mixing angle 
9 1 3 is small with sin 2 26> 13 < 0.1 [13], leading to further suppression factors as outlined in 
section 2.1. The smallness of these two parameters allows to classify the parameter space 
with the help of the analytic formulae in the scheme below, which is useful for reasons of 
structure and clarity. Our results are, however, based on a full numerical calculation, which 
essentially confirm this classification, but also fails in some places, as we will see later in 
this work. 

• Leading parameters: 

For Amj! = and 9 13 = 0, there are no transitions in the u e — > appearance channel. The 
disappearance probability reduces to the two neutrino case P{y^ — > v^) = 1 — sin 2 26> 2 3 sin 2 A 
being controlled by #23, Am^, which we call leading parameters. These parameters have 
already been measured by atmospheric neutrino experiments [14, 15, 16] and will be deter- 
mined better by conventional long baseline experiments [17, 18, 19, 20]. A neutrino factory 
will allow precision measurements of 6*23 and Ara^ and the result will be limited mainly by 
systematical errors. The measurement of these parameters will be dominated by the unsup- 
pressed rates in the — > and — > disappearance channels, which does not rely on 
excellent charge identification capabilities for secondary muons 2 . The differential event rate 
distribution allows precise fits of the energy spectrum and the question is only how good 
#23 and Am^ can be measured for certain detector and neutrino factory parameters, i.e. 
what the optimal beam energy and baseline are in this context. It will also be interesting 
to see how the accuracy for the leading parameters is modified for baselines which optimize 
the sensitivity to sin 2 2#i 3 or CP violation. This will be addressed in section 4.1. 

• Sub-leading parameters: 

For 9is 7^ and Am| x = 0, the first term in the appearance probabilities eq. (3) or in 
matter eq. (5) becomes non-zero: P[y e — > v^) = sin 2 #23 sin 2 2# 13 sin 2 ((A — 1)A)/(A — l) 2 . 
The appearance channel depends then via A on the sign of Am 31 [2, 3, 21]. In addition to the 
leading parameters, the analysis depends at this level also on #13 and sgnAm^ , which we call 
sub-leading parameters. Whether it is possible to determine these sub-leading parameters 
depends crucially on the value of #13, and the sensitivity limit below which no effects from # 13 
can be measured will be studied in detail. Matter effects have so far not been measured and 
an experimental test of MSW-effects [7, 8, 9, 10] will be possible if #13 is large enough. The 
measurement of #13 and the search for matter effects [2, 3, 1, 12] are thus important topics 

2 Sec e.g. [3] for a discussion of this problem. 
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for the physics program of a neutrino factory, which do not depend on Am^ and 9 12 being in 
the LMA-MSW range. Measurements of rates in the appearance channel depend, however, 
crucially on the capability to identify the charges of the secondary muons very reliably [3], 
which requires improved detector technology. Sensitivity limits to 6> 13 and statistical errors 
are studied in section 4.2, where also baseline and beam energy optimization are discussed. 

• Sub-sub-leading parameters: 

Finally, for Am^ 7^ and #13 7^ 0, effects due to the small solar mass squared splitting 
are added and the remaining three parameters appear in the oscillation formulae: Am^, 
9 12 and the CP phase #cp- Measuring leptonic CP violation is an exciting possibility for 
neutrino factories. In order to obtain sufficient rates, this requires, however, that the LMA- 
MSW region is the correct solution to the solar neutrino problem. One can see immediately 
from eq. (3) that Am^ (i.e. a) and #13 are the crucial parameters with determine the 
absolute and relative strength of the CP-violating effects. The foreseeable experiments are 
often not too far from the sensitivity limit. The limits on those two parameters, which 
will be discussed in section 4.3, are therefore among the most important points of this 
study, providing a better understanding how to search for CP violation and how 5 CP can 
be extracted. We will also see in section 4.3 that 9 12 and Am^ can be measured only 
very poorly in neutrino factory experiments. To a good approximation only the product 
Amjj sin 2# 2 i can be determined, which can be expected already from the the appearance 
probability in eq. (3). Thus we expect and will see that Am^ and 9 12 are highly correlated. 
Results obtained from neutrino factories cannot compete with results expected from the 
long baseline reactor experiment KamLAND [22]. It is thus very important to study how 
external input on Am^ and 9 12 (e.g. from KamLAND) helps to improve the extraction of 
CP-violating effects at a neutrino factory. The related questions about the optimal baseline 
and beam energy will be discussed in detail. As already mentioned above, the effects from 
the CP phase 5 C p are, like the matter effects, best accessible in the appearance channels 
v e — > Vfj, and D e — > which rely strongly on excellent charge identification capabilities of 
the detector. 

Note that the extraction of both the sub-leading and sub-sub-leading parameters is essen- 
tially based on moderate event rates, which can be estimated from the oscillation proba- 
bility at a mean energy. Unlike the disappearance channels, there is thus essentially no 
information in the energy spectrum. Note also that the above classification in leading, 
sub-leading and sub-sub-leading effects is also useful from a practical point of view. The 
experimental abilities of a neutrino factory will realistically develop from an initial, low flux 
setup, which is able to measure the leading parameters only, to a high luminosity machine, 
capable of measuring sub-leading and probably sub-sub-leading effects. This will allow in a 
good approximation an extraction of parameters at the relevant level, where the remaining 
parameters of the level below act essentially as unknown variables. 

2.3 Event rate analysis of neutrino factory experiments 

Our numerical analysis (see ref. [3] for details) is based on simulated event rates, which 
are calculated using neutrino flux profiles, full three neutrino oscillation probabilities in 
matter, a realistic Earth matter profile, neutrino cross sections, detector mass, threshold 
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and energy resolution. Backgrounds and experimental uncertainties are not taken into 
account. The v e — , P e —, v^— and v^— fluxes produced by a symmetrically operated neutrino 
factory of given muon energy and polarization can easily be derived from the muon flux and 
the kinematics of muon decay. The transition probabilities are calculated numerically by 
propagating the neutrino states through a realistic Earth density profile. The resulting muon 
event rates are then computed by folding the neutrino flux, the transition probabilities and 
the cross sections with the energy resolution function of the detector. The assumed muon 
detection threshold is 4 GeV. The energy resolution of the detector is approximated by a 
Gaussian resolution function with a = 0.1E, which gives an energy resolution AE/E = 10%. 
This value resembles the energy resolution of MONOLITH [23] for z/ M charged current events. 

With this method, event rates are calculated for the v e — > z/ M , v e — > z/ M appearance and 
v \i v \xi — > disappearance transitions with 20 energy bins in each channel. We as- 
sume perfect muon charge separation, i.e. that the different appearance and disappearance 
channels can be well separated. This requires for the appearance channels excellent charge 
identification capabilities of the detector. Insufficient charge identification would reduce the 
physics reach of such an experiment as pointed out in [3] . The type of detector assumed for 
this study is a magnetized iron detector with a mass of 10 kt. With 2 • 10 20 muon decays per 
year we obtain a standard "luminosity" of m^t = 2 • 10 21 kt year. We consider, however, 
also variations of the parameter N^m^. Thus, different fluxes like an initial 10 19 muon 
decays per year or higher fluxes (or detector masses) are considered. The details of the sta- 
tistical analysis are explained in chapter 3. The analysis uses simultaneously the neutrino 
and anti-neutrino channels of both the disappearance and appearance rate vectors, leading 
in all cases to optimal results. The method works of course also in cases where the results 
are dominated by one channel. Note, that is the sum of muons of both polarities, which 
means that the parameter m kt normalizes the sum of the neutrino and the anti-neutrino 
channels. Further details of the simulations can be found in ref. [3]. 

3 Statistical methods 

There are a number of non-trivial aspects in neutrino factory studies which require suitable 
statistical methods. We describe therefore in this section in more detail the problems and 
the methods used in this work. We also compare our method with other previously used 
methods. 

The aim of the analysis is to obtain statistically reliable statements about possible measure- 
ments of Am^, Am| l5 #12, #13, #23 and Sqp- A numerical scan of the full six dimensional 
input parameter with a subsequent extraction of the input parameters is an enormous nu- 
merical task. Therefore, up to now most studies performed only two parameter fits to 
evaluate the physics reach of a neutrino factory [12, 3, 24, 21, 25, 26, 27, 28, 1]. Often two 
parameters of the full parameter set, which were thought to be most relevant, were selected 
and a fit to this two parameters was performed using a x 2 -method. All other parameters 
were at the same time fixed to their best fit (i.e. input) values. Typically Gauss or Poisson 



6 



X 2 -functions have been used, i.e. 



Gauss x 2 = (6) 



i=i 

6 



Poisson x 2 = ^('2[(x i )-a; i ] + 2x i log 7 ^- ) . (7) 
i=i ^ 



{Xi) 



where b denotes the number of bins, X{ is the number of events in bin i and (x^ denotes the 
expectation value for bin i. o"j is estimated by yjxl. This method is numerically straightfor- 
ward, but it is not easy to include uncertainties of the left out parameters and systematical 
errors. The Gaussian x 2 -function allows, however, an approximate treatment of such un- 
certainties by Gaussian error propagation: 

X 2 = J> - <*<» ViJ 1 - , (8) 



d(xj) d{xj) 2 
, d\ a d\ a 

Here \ a is one out of n (in our case n — 6) parameters and a a denotes the uncertainty of 
A a . This method was for example used in [29]. It is, however, important to keep in mind 
that such an approach relies on the validity of the linear approximation of the functional 
dependence of [xj) on A. The point is that Gaussian error propagation relies only on the 
linear terms of a Taylor expansion of (xj) in A. The above approach may thus be no longer 
justified as soon as the second or higher order terms become important. 

Using one of eqs. 6-8, the acceptance region M at confidence level a is given by: 

X 2 (A) < X 2 (A) VA, (9) 
A X 2 (A) := X 2 (A)-X 2 (A), 

M = {\:A x 2 (\)<Axl}. (10) 

Using the asymptotic properties of maximum likelihood estimators, is determined from 
the condition that the area under a ^-distribution from zero to equals a. The number 
of degrees of freedom for this distribution is given by the number of free parameters r in 
the fit. The error on the confidence level a is then given by the size of the r-dimensional 
manifold 3 M. This procedure leads to an error interval I\ t for each of the r fitted parameters 
by simply projecting M onto Aj. We define the infimum and supremum with respect to \, 
i.e. the lower (upper) bound on A^: 

inf(M) := inf({Ai : A G M}) , 

Ai 

sup(M) := sup({Ai : A G M}) . (11) 

Ai 



3 Note that M is not necessarily connected and it may have a quite complicated topological structure, 
especially if there are parameter degeneracies. 
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Figure 1: Illustration of the error determination for a single parameter (in this example 
sin 2 26 ) i 3 ). In the upper row sin 2 29i 3 is fitted against the remaining 5 other parameters (grey 
regions) and the individual error bars of our method for each parameter are overlayed. These 
error bars result from the projection of the grey regions onto sin 2 2^13. All these errors are 
then combined into a total error for sin 2 2#i3, given by the lowest and highest endpoints 
of the individual errors. The result is shown as a grey region in the lower row, where the 
individual errors are overlayed. 



I\ i is then given by: 

(12) 

This method leads for n = 2 and a Gaussian likelihood to the usual error ellipses. 

We use a method which we call the Generalized Covariance Matrix Method (gcmm), 
where x 2 is calculated according to eq. 7 with n = 6 parameters, namely A = (Am^, Am^, 
#12, #13, #23, #cp)- The computational effort required to calculate the full six-dimensional 
manifold M is enormous 4 . To reduce this effort we adopted a scheme, which takes all 2- 
parameter-correlations into account. We calculate therefore all two dimensional slices of 
M through A which are parallel to the coordinate planes in order to estimate the errors for 
all six parameters: 

Pij = {A : A fc = X k ;k ^i,j}, 

= {A e : A% 2 (A) < Axi} ■ (13) 

These slices can be seen in the upper row of fig. 1, where the method is graphically explained. 
There are five such slices for one parameter and the largest extension of the projection 
of these five slices down to the parameter of interest is the corresponding error of this 

4 Even with our approximation below this analysis required to compute > 10 8 event rate spectra for each 
channel, which amounts roughly to l,000h CPU time on a 1 GHz Pentium III processor 



/a, = 



inf(M),sup(M) 
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parameter. The error interval I\ i of a parameter Aj is thus given by 



U 



mi(S ik ),sup(S ik ) 



(14) 



which corresponds to the grey band in the lower row of fig. 1. We use for Ax& the value 
9.2 which corresponds to 99% C.L. for two degrees of freedom, and to 83.8% C.L. for 
six degrees of freedom. The big advantage of the gcmm method is that no choice about 
"most influential" parameters has to be made. Instead, the method automatically finds the 
most influential parameters, like in the example of fig. 1, where the biggest error for # 13 
surprisingly comes from Am^. Previous neutrino factory studies used so far only one of 
the Sij to estimate the accuracy to which parameters can be extracted. We will show in 
this work that this can be a poor approximation, since some of the parameters, like #13 and 
Am^, are strongly correlated in a previously overlooked way. We will see that the sensitivity 
to sin 2 2#i3 can be reduced in this way up to one order of magnitude. Another advantage 
of the GCMM method is, that it can cope with uncertainties of all parameters even if the 
parameter dependence is highly non-linear, as it is often the case in neutrino oscillation 
studies. The GCMM method is thus in a number of ways better than a Gaussian error 
propagation. The GCMM is, however, still only an approximation, and full six parameter 
fits will lead to corrections. We compared the GCMM method to three parameter fits and 
found that it works quite good, but the difference can in the worst case amount to a factor 
of two. 

A final aspect of the analysis concerns the inclusion of parameters from other experiments. 
This leads to an understanding of how the parameter improvements of other experiments 
affect the analysis. An important case is Am^, which can be measured by KamLAND [22] 
in the LMA-MSW case with a precision which is much better than what can be obtained at 
a neutrino factory, as will be discussed in detail in section 4.3. It is therefore important to 
perform an analysis, where such external information can be included and where its impact 
can be assessed. Usually we extract in our analysis parameters by finding their best fit 
value and by determining the errors with the fitting procedure to our "data" . In order to 
include external knowledge of some parameter (e.g. on Am^), we simply restrict the range 
of variation of the is parameter 5 in the fits to the error interval at confidence level a. 

For the CP-phase there will, however, exist no information in advance. Since Sep determines 
which correlations are important, the errors on specific parameters can significantly depend 
on the value of 5qp- In such cases, the analysis is done for all possible values of 5qp and 
the maximal error which appears is taken as the final error. We call this procedure u Sqp 
unknown" . It is equivalent to integrating out a nuisance parameter with uniform prior in a 
Bayesian analysis. 

5 In a strict frequentist sense this leads to the loss of coverage. It should be understood as a uniform 
prior in the Bayesian sense on this parameter. 
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4 Results 



The results obtained from the numerical study are presented following the outline given 
in sec. 2.2. First, the leading parameters 6*23 and Am^ are discussed in sec. 4.1. Then, 
measurements of 13 , matter effects and the sign of Am^ are studied in sec. 4.2. Finally, CP- 
violating effects and measurements of Am| x and 6 12 are presented in sec. 4.3. The reader can 
select certain topics without missing information which is necessary for the understanding 
since the contents of these sections widely do not depend on each other. 

4.1 Leading parameters 2 3 and Am^ 

First, we will discuss the statistical sensitivity to the leading parameters #23 and Am^. Neu- 
trino factory experiments will allow high precision measurements of these two parameters, 
which are at present only roughly determined. The information on these parameters domi- 
nantly stems from disappearance channel measurements. Thus, the dependence on the value 
of #13 is marginal, especially when #13 is small. It will be demonstrated that the influence 
from uncertainties on Am^ and #21 can however be substantial. This has the consequence 
that input from other experiments (like KamLAND) can be helpful to limit the errors of 
the affected measurements. 

Optimization of baseline L and muon energy E^ 

Fig. 2 shows the expected statistical error on the quantities Am^ and sin 2 26* 2 3 in depen- 
dence of the neutrino factory parameters, baseline L and muon energy at N^m^ = 
2 • 10 21 kt year. For beam energies above 20 GeV and baselines above 2000 km, the statis- 
tical error on Am| x takes values between roughly 8% and 16%. For not too small muon 
energies, the important constraint is the baseline limit. The above mentioned 2000 km define 
the limit where the accuracy loss relative to the best value is less then a factor of two. This 
baseline limit linearly depends on the inverse of Am 2 1? which was chosen to be 3.5 • 10~ 3 eV 2 
in the above case. For Am^ = 1.0 • 10 _3 eV 2 the limit shifts to approximately 5000 km 
above which the statistical error is less than 20%. For Am 31 = 6.5 • 10~ 3 eV 2 the limit is 
well below 1000 km and the statistical error does not exceed 14%. 

The statistical error in measurements of the mixing angle #23 (right plot of fig. 2) is, with 
minimally 6%, at the same level of the one on Am^. Also here, a baseline limit can 
be given under which the statistical error is larger than twice the best achievable value. 
For Am 2 ! = 3.5 • 10~ 3 eV 2 (1.0 • 10~ 3 eV 2 ) this limit is roughly 3000km (5000km). The 
figure indicates that there is a favored value of L/E^, which can easily be understood using 
analytic considerations: Information on #23 is dominantly extracted from the total rates 
observed in the disappearance channel Piy^ — > u^) which is is not significantly affected by 
matter effects. Assuming a discrete neutrino energy corresponding to the average neutrino 
energy in the disappearance channel ((E v ) = 16/21£' M ), the number of observed muon 
neutrino events is approximately given by N = (1 — sin 2 2#23 sin 2 A31 )<!>„ , where the factor 
<E>„ M includes flux, cross-sections and detector mass. For = 2 • 10 21 kt year, we find 

that $^ = 6.910 6 (L/km)" 2 (16/21 J B M /GeV 2 ) 3 . In the Gaussian limit, which is a good 
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Figure 2: Statistical (relative) error of Am^ (left plot) and sin 2 2# 23 (right plot) as function 
of the baseline L and the muon energy for Am^ = 3.5 • 10 _3 eV 2 , sin 2 2# 23 = 1 and 
N^m^t = 2 • 10 21 ktyear. Dark shading indicates the preferred regions. For <5(Am 31 ), the 
best sensitivity reached is <5 min = 7.7%. For 5{/S.m\ l ) it is 6.1%. The contour lines indicate 
factors of 2, 5 and 10 accuracy loss compared to the best values. The sub-leading parameters 
sin 2 26* 13 , Amjj and 5cp do not play a crucial role and are assumed as unknown (see sec. 3). 

approximation here, the relative statistical error on the quantity sin 2 2# 23 is given by 
A sin 2 2023 



sin 2 26 23 



1 y/l- sin 2 2 6 23 sin 2 A 31 ~ _ L 1/2 

sin 2 2fl 23 sin 2 A 31 E^ ' 1 ' 



Thus, the relative error on sin 2 2# 23 is approximately a function of L/E^ with a slight 

— 1/2 

modification from the factor E^ . For muon energies between 10 GeV and 50 GeV this 
modification maximally gives a correction of (50/10) _1 / 2 ~ 0.45 which in general reduces the 
error for higher muon energies. The corresponding plot in fig. 2 clearly shows the contour 
lines of constant L/E^. Using eq. 15 with the parameters corresponding to fig. 2 yields 
a relative error for sin 2 2# 23 at the level 1/1000. For the quantity # 23 this translates to a 
relative error at the level of some per cent, which is in good agreement with the results 
obtained in the full numerical analysis. 

Summarily, the leading parameters do not give very strong recommendations for the se- 
lection of the baseline and muon energy. Beam energies between 30 GeV and 50 GeV are 
recommended. To achieve a sufficiently developed oscillation, the baseline must be large 
enough. For the central value of Am^ proposed by the Super-Kamiokande experiment, 
a baseline of 3000 km or more would be suitable. Should it turn out in the future, that 
Am 31 is at the lower edge of the presently favored region, one should be aware that this 
baseline limit shifts to higher values roughly inverse proportionally to Am|j. The exact 
values of the sub-leading parameters 6*i 3 and #cp do not play a role in this discussion since 
the results presented here depend only weakly on them. But Am 21 being correlated with 



11 



Am^ can have some impact on the errors made in measurements of Am^. This point will 
be discussed in detail later. 

Dependance on flux and detector mass iV M m^t 




Luminosity [2 N^m kt ] Luminosity [2 N^m kl ] 

Figure 3: Statistical error of 623 (left plot) and Am^ (right plot) as function of m^t, the 
number of stored muons per year times the mass of the neutrino detector in kilotons. The 
shaded bands indicate the range which is covered by variations of the sub-leadingparameters 
013 and Am§! (0 < sin 2 26 13 < 10 _1 , 10~ 5 < Am^ < KT 4 ). The calculation of 5(6 23 ) was 
performed with 50 GeV muon energy and a baseline of 7000 km (6823) respectively 3000 km 
(5Am 2 31 ). 

The influence of , the number of stored muons per year times the mass of the neutrino 

detector in kilotons, is shown in fig. 3. The decrease of the statistical error with increas- 
ing N^m^t roughly follows the prediction from Gaussian statistics (1/ y/N^ m^), which is 
indicated by the parallel lines in the background of the figure. The spread of the shaded 
band is generated by variations of the sub-leading parameters. Their influence increases 
with higher luminosities where the statistical errors are small. 

Correlation between Am^ and Am^ 

It is important to note that, for specific selections of L and E^, the statistical error on 
Amgj is dominated by its correlation with Am^. This fact is displayed in fig. 4, which 
shows the results of two parameter fits of Am^ with all other oscillation parameters. The 
plot shows that the dominating contribution to the statistical error of Am^ stems from the 
uncertainties of Am^ and 6u- This suggests that in such cases, external information on 
Anijj and #12 from other experiments can considerably improve the precision of measure- 
ments of Am| x . With fig. 5, this point is investigated in more detail: There, the resulting 
statistical error on Am^ is given as function of ^Am^), the error on Am| x as provided 
from external measurements. The shaded band indicates the expected values of ^Am^) 
from measurements performed at the KamLAND experiment [22]. The plots show that 
KamLAND input can improve the statistical error on Am^j up to a factor of three. This, 
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Am 2 2i Am 2 3 i 012 sin 2 20 13 e 32 <5cp 

Figure 4: Error bars of all two parameter fits which include Ara^. The total error (grey 
shaded region) is dominated by the contributions from the fits against Am^ and # 12 . Fit- 
ting Amjj only against 823 can result in a considerable underestimation of 5(Am^ 1 ). The 
calculation was performed with L = 8000 km, = 50 GeV, 9 12 = vr/4, sin 2 29 13 = 10" 3 , 
Amli = 10~ 4 eV 2 and iV^ m kt = 2 • 10 21 kt year. 

however, is only valid for specific baselines. At smaller baselines this correlation is not very 
pronounced. 




10° 10 1 10 2 10° 10 1 10 2 

(5Am 2 2 i (5Am 2 2 i 



Figure 5: Statistical error on Am 31 as function of ^Am^), the error on Am^ as provided 
from external measurements (like the KamLAND experiment). The shaded band indicates 
the expected values of S^Am^) from measurements performed at the KamLAND experiment 
[22]. No external input corresponds to the right edge of the plot (S(Aml 1 ) = 10 2 ). The 
calculation was done for Am^ = 10~ 4 eV 2 (left panel) respectively for Am^ = 10~ 5 eV 2 
(right panel) and 9 12 = vr/4. The different line styles indicate different values of sin 2 2# 13 

(io~ 4 , 10- 3 , 10~ 2 , lO^ 1 ). 
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4.2 Sub-leading parameters 6 13 and sgnAra 
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31 



Following the outline given in section 2.2, we continue this study with a discussion of 
the parameter # 13 . Information on 6* 13 is mainly obtained from the appearance channel 



probabilities P( 



and P(v P 



which contain terms proportional to sin2# 13 



and sin 2 2^13. For not too small values of sin 2 2#i3, the results presented in this section are 
roughly independent from the sub-sub-leading parameters Am^, #12 and <5cp. In this case, 
a neutrino factory experiment will be able to measure #13 with some precision. However, for 
small values of sin 2 2#i 3 close to the sensitivity reach, the parameters Am^, Qyi and <5cp can 
have considerable impact on measurements of 6*i 3 . This is an example where the step-wise 
analysis of such experiments is no longer valid. 



Optimization of baseline L and muon energy 
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Figure 6: Statistical error of the quantity sin 

sin 
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26* 13 as function of the baseline L and the 
muon energy for Am 31 = 3.5 • 10 -3 eV 2 , sin 2 2# 23 = 1, = 2 ■ 10 21 kt year and three 
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Dark shading indicates the preferred regions. 



The 



parameters Am^ and #cp play a role mainly for small values of 8i 3 . Here, they are assumed 
as unknown (see sec. 3). 
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Table 1: Statistical errors 5 made in measurements of sin 2 26*i 3 listed for several values 

The quantity 5, whose values are 
detailed description of 5 is given in the text. 



of Am 31 . The quantity 5, whose values are given in the table is the same as in fig. 6. A 



The statistical error of sin 2 2#i 3 in dependence of the neutrino factory parameters baseline 
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L and beam energy is displayed in fig. 6. The plotted quantity is 

5 := max{(2 min /z )" 1 ,W 2 o} , (16) 

where z = sin 2 2#i 3 and the subscripts correspond to the central value (0), the minimal 
value (min) and the maximal value (max) which are compatible with the simulated ex- 
perimental data. For values of 8 close to 1, 5 — 1 is approximately equal to the relative 
error. 5 = 10(100) indicates one (two) order(s) of magnitude uncertainty on the measured 
quantity. Since the error strongly depends on the value sin 2 2#i 3 , three plots for the values 
sin 2 26 l i3 = (10 _1 , 10^ 2 , 10~ 3 ) are given. For sin 2 2#i 3 = 10 -1 (left plot), the sensitivity is 
good and reaches values between 10% and 20% relative error at baselines between 3000 km 
and 8000 km and with beam energies above 20 GeV. For sin 2 2#i 3 = 10 -2 (middle plot), 
still errors down to 50% are reachable. Close to the sensitivity reach sin 2 2# 13 = 10" 3 (right 
plot), the situation is interesting. There, the information is still sufficient to determine 
sin 2 2# 13 with an error corresponding to roughly one order of magnitude. But now, close to 
the sensitivity limit, long baselines between 7000 km and 8000 km are strongly preferable. 
The reason for this is, that for small values of sin 2 2# 13 , the correlations of 9 13 with Am^ 
and 5 C p are important. These correlations tend to have less influence at large baselines. 
Since previous studies concerning the optimization of the baseline [12] did not include these 
correlations, they missed this point and came to differing recommendations for the baseline. 
Note, that close to the sensitivity limit, the result depends strongly on the precise value 
of the CP phase <5 C p- For the calculation which is presented here, we assume that 5 C p is 
completely unknown: The statistical errors are computed for several values of 5qp 6 and the 
maximal error of this set is considered as the final result. 

In table 1 the statistical errors made in measurements of sin 2 2# 13 are listed for several values 
of Ara^. The value of Am^ changes the obtained results significantly. Higher values of 
Am 31 increase the precision, whereas lower values make it more difficult to measure #i 3 . 

Correlation between 13 and the sub-sub-leading parameters Am^ and <5 C p 

In connection with measurements of <5cp it has been recognized that the correlation of the 
CP phase 5qp with the mixing angle 9 13 could possibly hide the effect of CP violation [29]. 
The implications of this correlation on the expected precision of measurements of of 0i 3 
have not been discussed. Already in the baseline discussion above, it was shown that in 
case of small values of 13 , correlations with 5 C p and Am^ can drastically influence the 
statistical errors expected in measurements of sin 2 2#i 3 . To illustrate this, fig. 7 shows the 
results of two parameter fits of sin 2 2#i 3 with all other oscillation parameters for two different 
values of the CP phase ^cp- The plots demonstrate that the dominating contribution to the 
statistical error of sin 2 29 13 can be caused by correlations with the parameters Ara^ or 5 CP . 
In the case 5qp = (left plot), the main error stems from the correlation with Ara^ but for 
5qp = n/2 (right plot), it is the CP phase 5qp which gives the dominating contribution to 
the total error. The correlation with Ara^ severely affects the sensitivity to sin 2 2#i 3 which 
was completely overlooked up to now. 

6 uniformly covering an interval of 2tt 
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Figure 7: Error bars of all two parameter fits which include sin 2 2#i3 for <5cp = (left 
plot) and Sqp = 7r/2 (right plot). The total error (grey shaded region) is dominated by the 
contributions from the fits against Am^ (left plot) or 5qp (right plot). Fitting sin 2 26> 13 only 
together with Ara^ can result in a considerable underestimation of the statistical error. The 
calculation was performed with L = 2100 km, = 50 GeV, 12 = tt/4, sin 2 29 13 = 10~ 2 , 
Am 2 x = 10~ 4 eV 2 and iV M m kt = 2 ■ 10 21 kt year. 



The influence of the sub-sub-leading parameters Ara^ and Sqp for the measurement of #13 
can qualitatively be understood using analytic considerations. For small values of #13 at 
the sensitivity limit, all four terms of equation 5 are equally important. The first, leading 
term has always a positive sign. The second term oc a sin Sep sin 2 2013 has a sign which 
depends on 5qp and whether neutrinos or anti-neutrinos are considered. Since we always 
use neutrinos and anti-neutrinos, this term can have both signs. The sign of the third term 
oc a cos^cp sin 2 2# 13 is determined by the value of Sqp- The sign of the fourth term oc a 2 
is always positive. If cos^cp < an increase (decrease) in 9 13 can be compensated by an 
increase (decrease) in Am^. In this case the main problem arises from the correlation of 6*13 
and Sqp, since all terms containing 5 CP are proportional to Am 2 x sin 2 2# 13 and this product 
can not be determined very well. Also the event rates are considerably smaller in this case 
which simply leads to a loss of statistical significance especially close to the sensitivity limit. 
If cos#cp > an increase (decrease) in 6>i 3 can be compensated by a decrease (increase) 
in Am^, therefore the product Am 2 x sin 2 26>i 3 can be determined quite accurately. In this 
case where the correlation of sin 2 26*13 with Am^ causes the problems, external information 
on Am^ (e.g. from KamLAND) would help to improve the results. 

It is also possible to understand, why large baselines around 7000 km are recommended in 
case of low values of #13 and large values of Am^: The relative size of the terms containing 
5qp decreases with growing baseline. Thus in principle larger baselines should perform 
better. For too long baselines, the l/L 2 -decrease of event rates starts to worsen the result 
again. The numerical calculation shows that baselines around 7000 km perform by far best. 
Shorter baselines can profit considerably from external input on Am 2 x if (and only if) 
cos5 C p > 0. 

It is of great importance to recognize that for small values of #13 and large values of Am 2 x 
the contributions to the total error of measurements of #13, originating from the sub-sub- 
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leading sector, are substantial and not only minor corrections. With decreasing Am^, 
these contributions get smaller and in the limit Am| x = there is no influence left from the 
parameters Am^, 9 12 and 5 C p- The same is true for small values of 9 12 (SMA solution). 



Sensitivity reach for measurements of 6 1S 




Figure 8: Sensitivity reach for measurements of sin 2 26*13. The area to the left of the lines 
indicates the parameter range where measurements are compatible with sin 2 26*13 — at 
99% C.L. . The calculation was performed with a beam energy of 50 GeV. The different line 
types are for different baseline as explained in the legend. For comparison to older studies, 
the left panel displays the result obtained from a two parameter fit of only sin 2 26?i3 and 
Am^. The right panel displays our new result with all parameters taken into account. 

We define the "sensitivity reach" for measurements of the mixing angle #13 as the max- 
imal value of sin 2 2O13 at which the experimental data is compatible with the hypothesis 
sin 2 2#i3 — at 99% C.L. . Measurements of nonzero values for #13 are only possible above 
the sensitivity reach. The sensitivity reach as function of Am^ is given by the right hand 
plot in fig. 8. The different line types indicate different baselines. For very small Am 2 ,!, 
values down to sin 2 26*13 — 10~ 4 are reachable and baselines between 1000 km and 7000 km 
perform very similar. The sensitivity limit varies only within a factor two. With increasing 
Ara^, the effects from the parameters Am^ and Sep become stronger and the sensitivity 
reach deteriorates by more than one order of magnitude (for short baselines). The left hand 
plot shows the results obtained by a two parameter fit to only #13 and Ara^. Since there, 
the correlation with Am^ and <5cp are not taken into account and <5cp is fixed to zero, the 
performance of an experiment is nearly independent of Am^. Comparing the two plots, 
it is obvious that with this method, for large values of Am^, the performance would be 
strongly overestimated, especially for baselines below 7000 km. This result is important 
since it suggests to use a longer baseline for measurements of 6*13 than usually proposed. 

The dependence of the #13 sensitivity reach on the parameter N^m^t is shown in fig. 9. 
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The scaling behavior is nearly Gaussian, i.e. it is roughly proportional to a/ . At 
high luminosities systematical errors, like insufficient charge identification capabilities of the 
detector, will most probably limit the sensitivity. 
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Figure 9: Sensitivity reach (see text for detailed definition) for measurements of sin 2 2# 13 
as function of m^, the number of stored muons per year times the mass of the neutrino 
detector in kilotons. The band indicates the range which is covered by variations of the sub- 
leadingparameters Am^ and Sqp (10~ 5 < Am^ < 10~ 4 , — tt < 5qp < 7r). The calculation 
was performed with a baseline of 7000 km and 50 GeV muon energy. The vertical dashed 
line indicates a typical neutrino factory experiment with a 10 kt iron detector. 



Sensitivity reach for the determination of sgn Am^ 

The ability to determine the sign of Am^ stems from the fact that the MSW-resonance 
occurs - depending on the sign of Am^ - either in the neutrino or in the anti-neutrino 
appearance channel. Therefore it is very important to be sensitive to the energy region 
around the MSW-resonance. In references [2, 3, 28] it was shown that, in the case this 
sensitivity is good, the ability to determine the sign of Am^ nearly coincides with the ability 
to determine #13. This holds in some approximation also when all parameters are taken into 
account. Hence, we do not study this point in detail but refer to fig. 9, which should give 
also a rough idea about the limit on # 13 above which a determination of sgnAm^ should 
be possible. Each term of eq. 5 is invariant under a simultaneous change from neutrino to 
anti- neutrino and of the sign of Am^. Therefore changing the sign of Am^ is equivalent 
to interchanging the role of neutrinos and anti-neutrinos. Assuming symmetric operation 
of the neutrino factory, the only difference that remains is due to the cross sections, which 
are for neutrinos twice as large as for anti-neutrinos. In case of a negative sign of Am^ one 
thus looses half of the statistics, which can be compensated by doubling the flux of u e . 



18 



4.3 Sub-sub-leading parameters 12 , Am^ and 6cp 

Whether it is possible to measure the parameters #i 2 , Am^i an d #cp with a neutrino factory 
long baseline experiment crucially depends on the values of #13, #i 2 and Ara^. To signif- 
icantly influence the measured event rates, it is absolutely necessary that the LMA-MSW 
solution is the correct description for solar neutrino oscillations. In comparison to the re- 
sults presented in the previous section, which do only depend on the value of #13, this is a 
severe additional condition which enables or disables the search for the CP-violating phase 
5 C p- Assuming now that the LMA-MSW solution is realized, it is important to notice that 
the sub-sub- leading effects are proportional to Am^, which in the LMA-range can still vary 
between roug hly 2 • lCT 5 eV 2 and 2 • 10~ 4 eV 2 (at 99% confidence level) [6]. This means, that 
the magnitude of CP-violating effects can also vary by one order of magnitude. 

In this section we will discuss the following points in detail: First, we will study the precision 
of measurements of the parameters 812 and Am^ in comparison to the KamLAND exper- 
iment. We will show, that the neutrino factory is mainly capable to measure the product 
Am| x sin 29 12 but has severe difficulties to constrain Am^ and 12 separately. KamLAND 
will provide much better accuracy for these parameters. Then, we will focus on the search 
for CP violation. We will discuss, for different neutrino factory fluxes and detector masses, 
the magnitude of the effects from the CP phase 5cp in dependence of the parameters # i3 and 
Am^. Which baselines and beam energies are best suitable for measurements of the CP 
phase is also studied in detail. Finally, some specific topics like the parameter degeneracy 
in the #i3-5cp parameter plane are discussed. 

Correlation between 12 and Am^ 

For Am 2 x < 10" 4 eV 2 , the appearance channel probabilities depend only on the product 
Arajj sin 2612 (see eq. 5). Thus a strong correlation between 612 and Am^ can be expected. 
The upper plot of fig. 10 shows the result of a two parameter fit of 9 2 \ an d Am^ to only the 
appearance rates. From the comparison to the dashed line, which represents the constant 
value Amj! sin2#i 2 = 10~ 4 eV 2 , it can clearly be seen that indeed the measurement is only 
sensitive to the product Am 21 sin26'i 2 . For large values of Am 21 , this correlation can be 
lifted by inclusion of disappearance rates (see lower plot of fig. 10). The sub- leading term 
in the disappearance probability (eq. 4) depends on the product Am^ cos 2 612, which helps 
to lift the degeneracy between Am 21 and 6*i 2 . If, however, Am| x is small, the a 2 -term in the 
disappearance probability looses in strength and it gets very difficult to obtain information 
on 9 12 and Am^ separately. 

Hence, the precision of measurements of 9 12 and Am^ is in general not very good. Particu- 
larly, the performance of a neutrino factory in such measurements is not comparable to the 
KamLAND experiment. Fig. 11 shows results of two parameter fits of 9u and Am 21 from the 
combined rates of the appearance and disappearance channels for several values of Am 2 x and 
9\2- The simulated experiment was performed with a luminosity of m^ t = 2- 10 22 kt year, 
which is by a factor ten higher than the standard value used in this work. It can be seen 
that despite the relatively high value of N^m^t, the results are not good. Particularly, the 
error ellipses grow drastically with decreasing 9 12 and Am 21 . 
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Figure 10: 3cr contour line of two parameter fits of 9 12 and Am^ to simulated event rates 
from the appearance channel only (upper plot) and appearance channel plus disappearance 
channel (lower plot). Constant values of the product Am^ sin2# 12 are given by the dashed 
lines. The calculation was performed with 50 GeV muon energy, 3000 km baseline, = 
2 • 10 21 ktyear, 9 12 = tt/4, Am 2 21 = 10" 4 eV 2 , sin 2 26 13 = 10~ 3 (upper plot) and sin 2 26 13 = 
10 _1 (lower plot). 

Correlation of <5 C p with 6 13 and possible degeneracies 

In section 4.2 it was demonstrated that the parameters #13 and <5cp can be strongly cor- 
related. This effect was visualized in the right plot of fig. 7. It is obviously important for 
studies of the CP-phase 5qp and it is discussed in detail in ref. [29]. Since our statisti- 
cal method takes into account all possible two parameter correlations, also this particular 
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Figure 11: Two parameter fits of 812 and Am^ to combined simulated event rates of the 
appearance channel and the disappearance channel. The calculation was performed with 
50 GeV muon energy, 3000km baseline, N^m^ t = 2 • 10 22 ktyear (!!) and sin 2 26*13 = 10 _1 . 

correlation is automatically included in our results. We do not further discuss this corre- 
lation but focus on another interesting problem which is also studied in the above referred 
work. There, it is demonstrated that considering the full parameter space of Sep, multiple 
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Figure 12: Fits to sin 2 2# 13 and Sqp with different values of the energy resolution of the 
detector (left plot) and x 2 -difference between the best fit points of the two degenerate 
solutions as function of the energy resolution (right plot) . The black square in the right 
plot denotes the standard energy resolution used in this work. Parameters: = 50 GeV, 
L = 3000 km, Am 2 x = 10~ 4 eV 2 . 

degenerate solutions appear in simultaneous fits of Sep and (9 13 . It is also stated that this 
degeneracy could possibly be lifted when different baselines or beam energies are studied 
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simultaneously. We find that the appearance of degenerate solutions crucially depends on 
the energy resolution of the detector. In fig. 12 the influence of the energy resolution on a 
fit of 5cp versus sin 2 29% 3 is shown. In the left hand plot the fit itself is depicted for the three 
energy resolutions 10%, 30% and 50%. The second degenerate solution in the upper part of 
the plot diminishes with increasing resolution and completely vanishes for the value 10%. 
To further illustrate this effect, the right hand plot shows the x 2 -difference of the best fit 
points of the two degenerate solutions as function of the energy resolution of the detector. 
The dashed horizontal lines in the plot indicate which energy resolution is needed to refuse 
the second degenerate solution at a certain confidence level. With an energy resolution of 
10%, like we use it throughout this study 7 , the second solution disappears with a confidence 
of more than 99%. We checked that this holds for all values of Sqp and #13. In ref. [29] the 
analysis was performed with only five energy bins. This corresponds to an energy resolution, 
which is too bad to allow a lift of the degeneracy. 

Optimization of baseline L and muon energy E^ 
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Figure 13: Results of fits of the CP phase <5cp as function of the baseline L and the 
muon energy for Am^ = 3.5 ■ 10~ 3 eV 2 , Am 2 21 = 10" 4 eV 2 , 6 23 = vr/4, 6 12 = vr/4, 
N^mu = 2 ■ 10 21 ktyear and three values of sin 2 2^13 (10 _1 , 10~ 2 , 10~ 3 ). Dark shading 
indicates the preferred regions. The quantity 5 plotted here is the percentage of the 5qp 
parameter space [— 7r/2,7r/2] which is compatible with the simulated experimental data at 
the 3a confidence level. The contour lines correspond to 5 = 50% and 5 = 95%. In the 
white shaded region no information on the CP phase can be obtained. 

In fig. 13 the regions in L-Z^-plane are shown where <5cp can be measured. We find that 
in general higher energies and baselines around 3000 km are preferred which is in good 
agreement with the results obtained in [12, 29] but is in clear contradiction to the results 
obtained in [30, 31]. The reason for this disagreement seems to lie in the different use of 
statistics. In the referred work no fits are performed but Gaussian error propagation is used 
(see sec. 3). In table 2 the optimal regions in the L-i? M -plane are given for different values 

7 This particular value seems reasonable since it is very close to the value which is quoted by the MONO- 
LITH collaboration [23]. Details to the definition of energy resolution can be found in section 2.3. 
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of sin 2 2#i3 . We also checked the influence of the value of Am^ on the sensitivity and the 
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Table 2: Optimal choice of the baseline L and the beam energy for measurements of 
the CP phase <5cp- The data of the first three lines are taken from figs. 13. The given 
ranges identify the regions where the statistical error 5(5cp) is not more than by a factor 
1.1 larger than the optimal value <5(<5cp) m in- In figs. 13 these regions are indicated by the 
darkest shading. 

optimal regions in the L-i^-plane (see table 3 ). The sensitivity gets better (worse) as 
Am 2 ! increases (decreases). The optimal baseline seems to be quite stable against changes 
of Am^. The optimal energy however gets smaller as Am^ gets smaller. 



Am 2 x 


sin 2 29 13 


fraction 


6.0- IO- 3 eV 2 
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22% 
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[ 28% 


3.5- IO" 3 eV 2 


io- 1 


26% 


IO' 2 


32% 


IO" 3 


40% 


1.0- IO" 3 eV 2 


io- 1 


100% 


io- 2 


100% 


IO- 3 


100% 



Table 3: Fraction of the parameter space of <5cp which is covered by the 3a acceptance 
region. 100% means that no information on the CP phase can be extracted from the 
experimental data. The values are based on Am 2 x = 10 _4 eV 2 , 9 2 3 = 7r/4, #i 2 = 7r/4 and 
m kt = 2 • 10 21 kt year. 



Sensitivity reach for Sep 

The two parameters Am^ and #13 control the size of all CP-effects. A very important 
question thus is, down to which values of these two parameters the CP phase can be deter- 
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mined. This sensitivity reach is shown in fig. 14. The three grey bands were obtained with 



10 



,-4 




flux [Npinn] 



| | 2 10 21 

□ 2 1022 
| 2 10 23 



10 



6 



10 



,-6 



10 



-5 



10 




,-3 



10 



,-2 



10 



,-1 



Figure 14: Sensitivity limit of measurements of the CP-phase 5qp in the O^-Am^ para- 
meter plane for three different values of N^m kt . The lower edge of each band represents 
the sensitivity limit under which no information on the CP-phase can be extracted from 
the experimental data. The upper edges of each band indicates 50% statistical error (at 
99% C.L.) for measurements of Sqp- This result was obtained with a baseline of 3000 km. 
The light grey shaded region in the background indicates the values of Am^ compatible 
with the LMA-MSW solution and the grey horizontal line indicates its best fit value. 

luminosities A^m kt of 2 • 10 21 ktyear (dark grey), 2 • 10 22 ktyear (grey) and 2 • 10 23 ktyear 
(light grey). The upper edge of each band represents 50% error, which is the contour line 
where half of the parameter space of <5cp can be ruled out by the experiment. The lower 
edge represents the sensitivity limit under which no information on 5cp can be extracted 
from the experimental data. The light grey shaded area in the background indicates the 
presently allowed LMA-MSW solution and the grey horizontal line its best fit value[6]. Note 
that the so called "initial stage" option 2 • 10 20 ktyear is not included in the plot since it 
does not provide sufficient event rates to make the CP-effects accessible. With the standard 
luminosity assumed in this work (2 • 10 21 kt year), CP-effects are only accessible for values of 
Araa at the upper edge of the LMA-solution. If Am| x is at the lower edge, the luminosity 
must be increased by a factor of more than 10. In this case, however, systematic errors 
could become the dominant source of uncertainties. 

5 Conclusion 

Previous neutrino factory studies were based on simplified treatments, which do not include 
all parameter correlations. We show that these correlations have considerable influence 
on the physics potential of a neutrino factory, especially for measurements close to the 
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sensitivity limit of #13 and for measurements of CP violation. We developed for our analysis 
improved statistical methods which are based on the computation of all two parameter 
correlations and which automatically include all parameter uncertainties. We used this 
method to look for all possible correlations and to refine in this way the understanding of 
the potential of neutrino factories. Using this method we found indeed correlations which 
were previously ignored or overlooked, leading to corrections for earlier results. The results 
which are summarized below were obtained with a default neutrino flux of 2 • 10 20 decaying 
muons per year directed to a 10 kt magnetized iron detector. Systematical, experimental 
and background limitations are not included in our study and the given errors and limits 
are thus of pure statistical nature. Unless mentioned differently the central value of the 
parameter Am^ was chosen to be 3.5 • 1(U 3 eV 2 . 

We found that the minimal errors which can be achieved for the leading parameters #23 and 
Airijj are roughly 8% and 6%, respectively. The precise length of the baseline is not crucial 
for these measurements, as long as it is larger than a certain minimal value. This lower limit 
depends strongly on Am^. For Am 2 x = 3.5 ■ lCT 3 eV 2 a baseline of 3000 km is enough to 
obtain the precision mentioned above. Lower values of Am| x required larger baselines. For 
Auigj = 1.0 • 10~ 3 eV 2 the baseline should, for example, exceed 5000 km to obtain optimal 
results and the precision which can be achieved is only 20%. Am 2 x values above 3.5- 10~ 3 eV 2 
allow on the other hand correspondingly shorter baselines. High beam energies are in general 
better for measurements of #23 and Ara^, but the energy dependence becomes rather weak 
between 30 GeV and 50 GeV. Correlations with the sub-leading parameters Am| x and 9n 
become for large Am^ important and contribute significantly to the total errors of #23 and 
Am^. A potential measurement of Ara^ in the LMA-MSW regime by KamLAND will 
therefore improve the errors on Am^ by up to a factor of three. 

For measurements of the mixing angle #13, we also found that beam energies between 30 GeV 
and 50 GeV perform approximately equally well, but two cases should be distinguished for 
the baseline. For 9i 3 somewhat below the current upper bound baselines between 3000 km 
and 8000 km are recommended. If, however, 9 13 is very small and therefore close to the 
sensitivity limit, we find a preferred baselines between 7000 km and 8000 km. The reason 
behind this is that there exists a quiet strong correlation between 9 13 and the sub-sub- 
leading parameters Am^, 9 12 and 5qp- This correlation is important for large values of 
Am 2 ! but it looses in strength for larger baselines. Our improved statistical analysis gives 
sensitivity limits for measurements of 9\ 3 which are up to one order of magnitude worse than 
previous results. If the LMA-MSW region is confirmed, then KamLAND's measurement of 
Am^ will improve the determination of #13 at a neutrino factory for cos 5qp > 0. We found 
that the sensitivity limit for 13 lies between sin 2 29i 3 = 10~ 4 and sin 2 29 13 = 10~ 3 , depending 
on the value of Am^ and on the baseline. A 10% statistical error on 9 n is expected if #13 
is close to the CHOOZ bound. A 50% error is still achievable for sin 2 29\ 3 = 10~ 2 . For a 
discussion of the flux dependence we refer to the corresponding sections of this work. We 
did not discuss in detail measurements of the sign of Am^ and the prove of the MSW-effect 
since these points are closely related to measurements of #13. One can thus roughly identify 
the sensitivity limits for measurements of 9 n with the limits down to which MSW-effects 
and the sign of Am^ are measurable. 

Finally, we focused on the parameters Am^, 612 and 5qp- In the LMA-MSW case a neutrino 
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factory can also measure Am^ and #12, but we have shown that a neutrino factory can not 
compete with KamLAND. To improve this situation, fluxes at least a factor 100 larger than 
usually discussed in the context of neutrino factories are necessary. The reason for this is 
that the appearance probabilities only depend on the product Am^ sin2#i 2 , which makes it 
difficult to obtain separately information on Am^ and #12. Concerning measurements of the 
CP-phase <5cp, we obtained the following results: Baselines between 2800 km and 3500 km 
are good choices. Lower baselines are not recommended, not only because of backgrounds 
which would spoil the signal, but also from a statistical point of view. Muon energies 
between 30 GeV and 50 GeV are again preferred for such measurements. We find, however, 
that N^m^t = 2 • 10 21 ktyear is lower bound for a CP violation measurement, which is at 
this limit only possible for very large values of Am^. To cover the full LMA-MSW region, 
a luminosity of at least m^ t = 2 • 10 22 kt year is needed. Moreover, we found that the 
degeneracy in the 13 -5qp parameter plane, which was recently pointed out in the literature 
is not present in our analysis. Such a degeneracy shows only up when the energy resolution 
of the detector is reduced. 

We would like to stress again, that we perform a statistical analysis and that systematic 
errors and backgrounds are not included. The adopted statistical method resembles to a 
large extent the results which would be obtained by full six parameter fits. Nevertheless, the 
errors are in some cases still somewhat underestimated and we found that full six parameter 
fits can have up to a factor two larger errors. 

Acknowledgments: We wish to thank E. Akhmedov and T. Ohlsson for discussions and 
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